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Abstract 

We consider a super symmetric (SUSY) scenario including right-handed neutrinos, 
one of whose scalar superpartners is the lightest SUSY particle (LSP). The distinguish- 
ing feature in the collider signal of SUSY in such a case is not missing energy but a 
pair of charged tracks corresponding to the next-to-lightest SUSY particle, when it is, 
as in the case considered, a stau. Following up on our recent work on neutralino re- 
construction in such cases, we explore the possibility of reconstructing charginos, too, 
through a study of transverse mass distributions in specified final states. The various 
steps of isolating the transverse momenta of neutrinos relevant for this are outlined, 
and regions of the parameter space where our procedure works are identified. 
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1 Introduction 



If supersymmetry (SUSY) [1, 2], broken at the TeV scale, has to validate itself as the next 
step in physics beyond the standard model (SM), then it is likely to be discovered at the 
Large Hadron Collider (LHC), with the superpartners of the SM particles being identified. 
It is therefore of great importance to make a thorough inventory of collider signals answering 
to various SUSY scenarios. These can serve not only to unveil the general character of the 
scenario but also to yield a wealth of information as possible about the specific properties of 
the new particles. 

One attractive feature of supersymmetric theories, in the i?-parity (defined as R = 
(— ) 3B+L+25 ') conserving form, is that the superparticles are always produced in pairs and 
each interaction vertex must involve an even number of superparticles (having R = -1). 
Hence the lightest SUSY particle (LSP) is stable and all SUSY cascades at the collider ex- 
periments should end up with the pair production of the LSP. As a bonus, it provides a viable 
candidate for dark matter (DM) if the LSP is electrically neutral and only weakly interact- 
ing. Since the LSP, due to such a character, escapes the detector without being detected, a 
prototype signature of i?-conserving SUSY is energetic jets and/or leptons associated with 
large missing transverse energy (fir)- 

However, there can be several SUSY models [3] which include a quasi-stable next-to- 
lightest supersymmetric particle (NLSP). This can happen in cases where the NLSP is nearly 
degenerate in mass with the LSP, or if its coupling to the LSP is too small. Hence the decay 
width of the NLSP into the LSP is highly suppressed. Consequently the NLSP becomes 
stable on the detector scale, its lifetime being long enough to escape the detector without 
decaying inside it. Thus the NLSP behaves like a stable particle within the detector. The 
resulting collider signals change drastically, especially if the NLSP is a charged particle. The 
quintessential SUSY signals then is not fir but two hard charged tracks of massive stable 
particles which appear as far as up to the muon chamber. This opens up a whole set of new 
possibilities for collider studies, including reconstruction of the sparticle masses, something 
that is relatively more difficult in the presence of fir- 

The scenario we have considered here, as an illustration of such a quasi-stable NLSP, is 
the minimal supersymmetric standard model (MSSM) augmented with a right-chiral Dirac- 
type neutrino superfields for each generation. This is consistent with the existing evidence 
[4] of neutrino masses and mixing, although no explanation for the smallness of neutrino 
masses is offered. It is possible in such a case to have an LSP is dominated by this right- 
chiral sneutrino state (v R ) together with a charged particle as the NLSP. We specifically 
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consider a situation with a stau (f) NLSP 1 . Such a scenario can easily be motivated [7] by 
assuming that the MSSM is embedded in a high scale framework of SUSY breaking. As we 
shall see, this can happen in minimal supergravity (mSUGRA) [8] where the masses evolve 
from "universal" scalar (m ) and gaugino (Mi/ 2 ) mass parameters at a high scale. The only 
extension here is the right-chiral neutrino superfield (in fact, three of them) whose scalar 
component derives its soft mass from the same m . The existence of such a quasistable 
charged slepton can be well in agreement with the observed abundance of light elements as 
predicted by the big-bang nucleosynthesis (BBN) [9], provided, its mass is below a TeV. Since 
the right-chiral sneutrino has no gauge couplling but only interactions proportional to the 
neutrino Yukawa coupling, the strength of which is too feeble to be seen in dark matter search 
experiments, such an LSP is consistent with all direct searches carried out so far. Moreover, 
it has been shown that such a spectrum is consistent with all low energy constraints [10], 
and the contribution to the relic density of the universe can be compatible with the limits 
set by the Wilkinson Microwave Anisotropy Probe (WMAP) [11] with appropriate values of 
the relevant parameters [12]. 

In this work, we have concentrated on the mass reconstruction of the lightest chargino 
in a f-NLSP, i>i-LSP scenario. This is a follow-up of our earlier work [13] on neutrlino 
reconstruction under similar circumstances. We have shown that it is possible to determine 
the mass of the lightest chargino, produced in the cascade decay of squarks or gluinos, from 
the sharp drop noticed in the transverse mass [14] distribution of the chargino decay products. 
More precisely, we show a way of disentangling the transverse mass of the system consisting 
of a f-track and the associated neutrino from chargino decay. We suggest a method for 
extracting the transverse momentum (p T ) of the neutrino. Though a sizable statistics is 
required for this purpose, and one may have to wait for considerable accumulated luminosity 
after the discovery of the LHC, still this is a rather spectacular prospect. We have successfully 
applied the criteria, developed in our earlier work, for separating the signal from standard 
model backgrounds. Ways of suppressing SUSY processes that are likely to contaminate the 
transverse mass distributions are also suggested. 

It should be mentioned that neither the signal we have studied nor the prescribed recon- 
struction technique is limited just to scenarios with right-sneutrino LSP. It can be applied 
successfully to all cases [15, 16, 17] where the NLSP is a charged scalar with quasi-stable 
character, provided that it decays outside the detector, leaving behind a charged track in 

1 It should be remembered that the above possibility is not unique. One may as well have a spectrum in 
which a third family squark is the NLSP [5, 6]. 
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the muon chamber. 

The paper is organized as follows: in Sec. 2, we motivate the scenario under investigation 
and present a brief review of the mass reconstruction for neutralinos as done in the earlier 
paper and used in this work. There we also summarise our choice of benchmark points. The 
signal under study and the reconstruction strategy for determining the chargino mass as well 
as the possible sources of background, both from the SM and within the model itself, and 
their possible discrimination, are discussed in Sec. 3. We summarise and conclude in Sec. 4. 

2 The overall scenario, mass reconstruction and repre- 
sentative benchmark points 

2.1 Scenarios with v R LSP and f NLSP 

As has been already stated, the most simple-minded extension of the MSSM [18], accommo- 
dating neutrino masses, is the addition of one right-handed neutrino superfield per family. 
In this situation the neutrinos have Dirac masses induced by very small Yukawa couplings. 
The superpotential of such an extended MSSM becomes (suppressing family indices), 

Wmssm = ViLH d E c + y d QH d D c + y u QH u U c + fiH d H u + y v H u Lv c R (1) 

where H d and H u , respectively, are the Higgs superfields that give masses to the T 3 = —1/2 
and T 3 = +1/2 fermions, and y's are the strengths of Yukawa interaction. L and Q are the 
left-handed lepton and quark superfields respectively, whereas E c , D c and U c , in that order, 
are the right handed gauge singlet charged lepton, down-type and up-type quark superfields. 
/i is the Higgsino mass parameter. 

It is a common practice to attempt reduction of free parameters in the theory, by assuming 
a high-scale framework of SUSY breaking. The most commonly adopted scheme is based on 
N = 1 minimal supergravity (mSUGRA). There SUSY breaking in the hidden sector at high 
scale is manifested in universal soft masses for scalars (m ) and gauginos (Mi/ 2 ), together 
with the trilinear (A) and bilinear (B) SUSY breaking parameters in the scalar sector. The 
bilinear parameter is determined by the electroweak symmetry breaking (EWSB) conditions. 
All the scalar and gaugino masses at low energy obtained by renormalization group evolution 
(RGE) of the universal mass parameters m Q and My 2 from high-scale values [19]. Thus one 
generates all the squark, slepton, and gaugino masses as well as all the mass parameters in 
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the Higgs sector. The Higgsino mass parameter \i (up to a sign), too is determined from 
EWSB conditions. All one has to do in this scheme is to specify the high scale (m , My 2 , A , 
together with sign(fx) and tan/3 = (H u ) / (Ha)) where, tan/5 is the ratio of the vacuum 
expectation values of the two Higgs doublets that give masses to the up-and down-type 
quarks respectively. 

The neutrino masses are typically given by, 

m v = y v (h° u ) = y v v sin/3 (2) 



The small Dirac masses of the neutrinos imply that the neutrino Yukawa couplings (y u ) are 
quite small 10~ 13 ). 

With the inclusion of the right-chiral neutrino superfields as a minimal extension, it 
makes sense to assume that the masses of their scalar components, too, originate in the 
same parameter m . The evolution of all other parameters practically remain the same in 
this scenario as in the MSSM, while the right-chiral sneutrino mass parameter evolves at the 
one-loop level [20] as 

_ 2 2 A * (3) 
dt ~ 16t* Vv u ■ [ ) 

where A v is obtained by the running of the trilinear soft SUSY breaking term A and is 
responsible for left-right mixing in the sneutrino mass matrix. 

It follows from above that the value of M„ R remains practically frozen at mo, thanks to 
the extremely small Yukawa couplings, whereas the other sfermion masses are enhanced at 
the electroweak scale. Thus, for a wide range of values of the gaugino masses, one naturally 
ends up with a sneutrino LSP (z>i), dominated by the right-chiral state. This is because the 
mixing angle is controlled by the neutrino Yukawa couplings: 

v\ = — u L sm9 + v R cos9 (4) 

where the mixing angle is given by, 

, OQ 2y v vsm P\ cot Pfi- A v \ 

tan2# = k k (5) 

Of the three charged sleptons, the amount of left-right mixing is always the largest in 
the third, and hence the lighter stau (fx) often turns out to be the NLSP in such a scenario. 
There are regions in the parameter space the three lighter sneutrino states corresponding 
to the three flavours act virtually as co-LSP's. It is, however, sufficient for illustrating our 
points to consider the lighter sneutrino mass eigenstate of the third family, as long as the 
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state (fi) is the lightest among the charged sleptons. Thus the addition of a right-handed 
sneutrino superfield, for each family, which is perhaps the most minimal input to explain 
neutrino masses and mixing, can eminently turn a mSUGRA theory into one with a stau 
NLSP and a sneutrino LSP. It is should be emphasized that the physical LSP state can 
have (a) Yukawa couplings proportional to the neutrino mass, and (b) gauge coupling with 
the small left-chiral admixture in it, driven by left-right mixing which is again proportional 
to y v . Thus the decay of any particle (particularly the NLSP) into the LSP will always 
be a very slow process, not taking place within the detector. Under such circumstances, 
the quintessential SUSY signal is not anymore but a pair of charged tracks left by the 
quasi-stable NLSP. 

2.2 Neutralino reconstructed 

In an earlier study [13], we suggested a reconstruction technique for at least one of the two 
lightest neutralinos, in the i>r LSP and f NLSP scenario. The signal studied there was: 
2tj + 2f '(charged — track) + Pt+ X. Here Tj denotes a jet out of one-prong decays of the tau, 
and all accompanying hard jets arising from cascades are included in X. The kinematic cuts 
imposed by us, such as p t r f ack >100 GeV and J2 \pt\ >1 TeV, (pr = transverse momentum, 
and YJ \pr\ is the scalar sum of all visible tarnsverse momenta) reduced the backgrounds 
considerably. 

Since the signal we investigated involves two taus in the final state, and hadronic decays 
of tau was considered, tau-jet identification and tau reconstruction were two important 
components of the procedure. For this a method suggested in [21] was used, which involved 
solving the following equation event-by-event: 

i>T = (— - 1) Phi + (— - 1) Ph2 (6) 

x Thl x Th2 

where x Thi (i = 1,2) is the fraction of the tau energy carried by each product jet collinear with 

— # 

the parent tau, when it is boosted. j>T is the vector sum of the transverse components of the 
3-momenta of the two product neutrinos produced in hadronic decay of each tau. Clearly 
this method is applicable when there is no other invisible particle in the final state. This 
was ensured in the best possible manner by vetoing any isolated lepton in the final state, 
thus getting rid of additional neutrinos from W-decays. 

— * — * 

The reconstruction of fa is undoubtedly very crucial here. j>T is reconstructed as the 
negative of the total visible p~j> which receives the contributions from isolated leptons/ slep- 
tons, jets and unclustered components. The last among these includes all particles (elec- 
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tron/photon/ muon/stau) with 0.5 < E T < lOGeV and \r]\ < 5 (for muon or muon-like 
tracks, \q\ < 2.5), or hadrons with 0.5 < E T < 20GeV and \rj\ < 5, which do not contribute 
to a jet and constitute 'hits' in the detector [22]. In order to simulate the finite resolu- 
tion of detectors, the energies/transverse momenta of all particles were smeared following 
prescriptions detailed in [13]. 

Once the tau 4-momenta are obtained, they are combined with the staus to find the 
stau-tau invariant mass distribution. This requires the knowledge of stau mass 2 as well as 
the choice of the correct pair. The correct pairs are obtained by using a seed mass for f 
which was taken to be 100 GeV, satisfying the criterion \M?f rl - M'^ ir2 \ < 50 GeV. The 
actual stau mass is then extracted by demanding that the invariant masses of the two fr 
pairs were equal, which yields an equation involving one unknown quantity, namely, m f : 

\jm 2 f + \p^\ 2 .E Tl - ^rr4 + \ P f 2 \ 2 .E T2 = p\.p~* Tl - Pf 2 .pr 2 (7) 

Having thus extracted mf on an event-by-event basis in the event generator, it was 
demonstrated that the distribution of this mass value has a peak at the actual mf. We have 
used this peak value in reconstructing the neutralino from the invariant mass distribution of 
the fr pair. In some regions of the parameter space, it is possible to thus reconstruct only 
Xi, as the production rate of m cascade decay of q (or g) as well as the decay branching 
ratio of X2 ~^ ^ r i s small. In some other regions, we have been able to reconstruct both 
of them. There are still other regions only the X2 peak shows up. This is because of the 
small mass splitting between Xi an d f , which softens the tau (jet) arising from its decay, 
preventing it from passing the requisite hardness cut. 

2.3 The choice of benchmark points 

The choice of benchmark points for this study is the same as in the case of neutralino 
reconstruction [13]. The mSUGRA parameter space is utilised for this purpose. A f NLSP 
and a i>\ LSP occur in those regions in which one would have had a f LSP in the absence of 
right-chiral neutrino superfields. We focus on both the regions where (a) m fl > + mw, 
and (b) the above inequality is not satisfied. In the first case, the dominant decay mode is 
the two-body decay of the NLSP, f\ — > uiW, and, in the second, the decay takes place via a 
virtual W. Decay into a charged Higgs is a subdominant channel for the lighter stau. The 
decay takes place outside the detector in all cases. At different benchmark points, however, 

2 From the muon chamber only the three-momenta of the charged track can be obtained. 
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the mass splittings between the fi and neutralinos/charginos are different. This in turn 
affects the kinematic characteristics of the final states under consideration. 

We have used the spectrum generator ISAJET 7.78 [23] for our study. In Table- 1 we list the 
six benchmark points used, both in terms of high-scale parameters and low-energy spectra. 
The justification of their choice and their representative character have been explained in 
reference [13]. 





BP-1 


BP-2 


BP-3 


BP-4 


BP-5 


BP-6 


mSUGRA 


m = 100 


m = 100 


m = 100 


m = 100 


m = 100 


m = 100 


input 


lTLi/2 = 600 


mi/2 = 500 


mi io = 400 


mi/2 = 350 


TTXi in = 325 


TTXi io = 325 




tan (3 = 30 


tan = 30 


tan = 30 


tan p = 30 


tan/? = 30 


tan p = 25 




418 


355 


292 


262 


247 


247 




246 


214 


183 


169 


162 


162 




408 


343 


279 


247 


232 


232 




395 


333 


270 


239 


224 


226 


m 0iR 


100 


100 


100 


100 


100 


100 


m fl 


189 


158 


127 


112 


106 


124 


fjif 2 


419 


359 


301 


273 


259 


255 


X-l 


248 


204 


161 


140 


129 


129 


m v o 


469 


386 


303 


261 


241 


240 




470 


387 


303 


262 


241 


241 


m- g 


1362 


1151 


937 


829 


774 


774 


m tl 


969 


816 


772 


582 


634 


543 


m t 2 


1179 


1008 


818 


750 


683 


709 




115 


114 


112 


111 


111 


111 



Table 1: Proposed benchmark points (BP) for the study of the stau-NLSP scenario in SUGRA with 
right- sneutrino LSP. The values of mo and M x / 2 are given in GeV. We have also set Aq = 100 GeV 
and sgn(fi) = + for benchmark points under study. 

It may be noted that the region of the mSUGRA parameter space where we have worked 
is consistent with all the experimental bounds [24], including both collider and low-energy 
constraints (such as the LEP and Tevatron constraints on the masses of Higgs, gluinos, 
charginos etc., as also those from b — > 57, correction to the p-parameter, (g^ — 2) and so 
on). In the next section we describe the procedure for the reconstruction of xt, f° r these 
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benchmark points. 



3 Reconstruction of the lighter chargino 

The final state of use for the reconstruction of the lighter chargino is 

• Tj + 2f(opposite — sign charged tracks) + fir + X 

where, Tj represents a jet which has been identified as a tau jet, the missing transverse energy 
is denoted by fi T and all other jets coming from cascade decays are included in X 3 . 

Simulation for the LHC has been done for the signal as well as backgrounds using PYTHIA 
(V6.4.16) [25]. The pp events has been studied with a centre-of-mass energy (E c . m ) = lA TeV 
at an integrated luminosity of 300 fb^ 1 . The numerical values of the electromagnetic and 
the strong coupling constant have been set at a~^(M z ) = 127.9 and a s (M z ) = 0.118 respec- 
tively [10]. The hard scattering process has been folded with CTEQ5L parton distribution 
function [26]. We have set the factorization and renormalization scales at the average mass 
of the particles produced in the parton level hard scattering process. In order to make 
our estimate conservative, the signal rates have not been multiplied by any K-factor [27]. 
The effects of initial and final state radiation as well as the finite detector resolution of the 
energies/momenta of the final state particles have been taken into account. 

3.1 Chargino reconstruction from transverse mass distribution 

Now we are all set to describe the main principle adopted by us for chargino (xt) recon- 
struction. For this purpose, we have looked for the processes in which xt ~ X1/X2 * s being 
produced in association with hard jets in cascade decays of squarks and gluinos. The xt 
subsequently decays into a f-v T pair, while the x\ ( or X2) decays into a f-r pair. Since the 
decay of xt involves an invisible particle (v T ), for which it is not possible to know all the 
four components of momenta, a transverse mass distribution, rather than invariant mass 
distribution, of r-v T pair will give us mass information of xt- m spite of the recent progress 
in measuring the masses of particles in semi-invisible decay mode (for example the % 2 vari- 
able introduced by [28] and its further implications [29]), we have focused on transverse mass 

3 In order to avoid the combinatorial backgrounds, we have considered events with two opposite sign f's 
only. 
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variable (tot) because of the fact that the only invisible particle present in the final state is 
the neutrino, which is massless. 

The procedure, however, still remains problematic, because the r on the other side (arising 
from neutralino decay) also produces a neutrino in the final state, which contributes to fix- 
In order to correctly reconstruct the transverse mass of the f — v T pair from chargino decay, 
the contribution to Ifi? from the aforementioned neutrino must be subtracted. 

Keeping this in mind, we have prescribed a method for reconstruction of the transverse 
component of the neutrino 4- momenta (Pj) produced from the decay of xt > i n association 
with f. To describe it in short: 

We label the transverse momentum of the neutrino coming from chargino decay by Pj. 
Attention is focused on cases where the tau, produced from a neutralino, decays hadronically 
and the r-jet, out of a one-prong decay of tau, is identified following the prescription of 
[21]. We have assumed a true tau-jet identification efficiency to be 50%, while a non-tau 
jet rejection factor of 100 has been used [30, 31, 32] (The results for higher identification 
efficiency are also shown in Sec. 4.). We have also assumed that there is no invisible particle 
other than the two neutrinos mentioned above. We have attempted to ensure this by vetoing 
any event with isolated charged leptons. This only leaves out neutrinos from Z-decay and 
PF-decays into a tu t pair. The contamination of our signal from these are found to be rather 
modest. 

The transverse momenta of the neutrino (Pj 2 ), out of a tau decay, is first reconstructed 
in the collinear approximation, where the product neutrino and the jet are both assumed to 
move collinearly with the parent tau. In this approximation, one can write 

P Tj = XP T (8) 

Following the decay Xi ( or ! xV)^ t^t t we have then combined the identified tau-jet with 
the oppositely charged stau (track), thus forming the invariant mass 

m 2 x , = {P f + P T ) 2 = {P f + P T Jx) 2 (2 = 1,2) (9) 

This pairing requires the charge information of the tau-induced jet. We have assumed 
that, for a true tau-jet, the charge identification efficiency is 100%, while to a non-tau jet we 
have randomly assigned positive and negative charge, each with 50% weight. One can solve 
this equation for x (neglecting the tau-jet invariant mass 4 ), to obtain 

4 This approximation is not valid for small x, say x < 0.1. However, the jet out of a tau decay almost 
always carries a larger fraction of r-cnergy, thus justifying the approximation. 
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2P f .P T , x 

x » r T3 (10) 
m o — w~ 

This requires the information of m x o (or, m x o) and as well, which we have used from 
our earlier work for the respective benchmark points. Once x is known we have, 

pT = pr _ pr = pT (n) 

Hence, the transverse component of the neutrino, out of xf-decay can be extracted from 
the knowledge of fix of that particular event 5 . This is given by, 

?l =$r- ?l (12) 

Finally, from the end point of the transverse mass distribution of the f-v T pair the value of 
m x ± can be obtained. However, one should keep in mind that, both Xi or X2 can decay into 
a tt pair. Therefore, it is necessary to specify some criterion to separate whether a given fr 
pair has originated from a x\ or X2 which we have discussed in the next subsection. 



3.2 Distinguishing between decay products of xi an d X2 

In order to identify the origin of a given opposite sign f-r pair, the first information that is 
to be extracted from data is which benchmark region one is in. We have assumed of gaugino 
mass universality for this process, for the sake of simplicity. 

— * 

If one looks at the effective mass (defined by M e ff = J2\pt\ + fir) distribution of the 
final state, then the peak of the distribution gives one an idea of the masses of the strongly 
interacting superparticles which are the dominant products of the initial hard scattering 
process. This is seen from Figure 1. Once the order of magnitude of the gluino mass is 
inferred from this distribution, one can use the universality of gaugino masses, which, in 
turn, indicates where m x o and m x o, masses of the two lightest neutralinos, are expected to 
lie. 

Next, for each event that we record, we look at a f and a r-jet of opposite signs. The 
invariant mass distribution of this ttj pair displays a peak whose location, although not 
precisely telling us about the parent neutralino, is still in the vicinity of the mass values. 
Thus, by observing these distributions (Figure 2) one often is able to tell whether it is a Xi 
or a x% once one simultaneously uses information obtained from the M e ff distribution. 

5 For details on the reconstruction of fir see [13]. 
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Effective Mass in GeV Effective Mass in GeV 

Figure 1: M e ff distribution (normalized to unity) of the final state under consideration, for all 
benchmark points. 

As has already been noted in [13] , the mass of either Xi or x% or both can be reconstructed 
in this scenario, depending on one's location in the parameter space. Once a peak in the tt 
invariant mass is located, the next step is to check whether \Mf_ T . — ?n x o| < 0.1.m x o, where 
m x o is either one (or the only one) of the two lightest neutralinos deemed reconstructible in 
the corresponding region. The mass of that neutralino is used in equation 9. If this equality 
is not satisfied for either neutralino or the only one reconstructed, then the event is not 
included in the analysis. 

3.3 SM backgrounds and cuts 

The final state we have considered, namely, Tj + 2f(opposite — sign charged track) + f, T + 
X, suffers from several SM background processes. This is because charged tracks in the 
muon chamber due to the presence of quasistable charged particle can be faked by muons. 
Such faking is particularly likely for ultra-relativistic particles, for which neither the time- 
delay measurement nor the degree of ionisation in the inner tracking chamber is a reliable 
discriminator. The dominant contributions come from the following subprocesses: 

1. tt: This is a potential background for any final state in the context of LHC, due to 
its large production cross-section. In this case tt — > bW + bW~, followed by various 
combinations of leptonic as well as hadronic decays of the b and the W, can produce 
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Invariant mass (stau-tau pair) in GeV Invariant mass (stau-tau pair) in GeV 




Invariant mass (stau-tau pair) in GeV Invariant mass (stau-tau pair) in GeV 



Figure 2: Mf Tj distribution for all the benchmark points. BP1, BP 2 and BP 3 show only the X\ 
peak. Both the Xi an d X2 P ea ks are visible for BP4 and BP5, while BP6 displays only the x% peak. 
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Muon pt in GeV Muon pt in GeV 



Figure 3: pr distribution (normalised to unity) of the harder muonlike track for the signal and 
the background, for all benchmark points. The vertical lines indicate the effects of a px-cut at 100 
GeV. 




Figure 4: £|pr| distribution (normalized to unity) for the signal and the background, for all bench- 
mark points. 
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opposite sign dimuons and jets. The jets may emanate from actual taus, but may 
as well be fake. One has an efficiency of 50% in the former case, and a mistagging 
probability of 1% in the latter. The it cross-section has been multiplied by a K-factor 
of 1.8 [33]. 

2. bb: This, too, has an overwhelmingly large event rate at the LHC The semileptonic 
decay of both the b's (b — > cfiu^) can give rise to a dimuon final state and any of the 
associated jets can be faked as tau-jet. Though the mistagging probability of a non- 
tau jet being identified as a tau-jet is small, the large cross-section of bb production 
warrants serious attention to this background. 

3. ZZ : In this case any one of the Z's can decays into a dimuon pair (Z — > fifi) while 
the other one can decays into rr pair where only one of the tau can be identified. The 
hadronic decay of Z and the subsequent misidentification of any of them as tau-jet is 
also possible. 

4. ZW : This SM process also contributes to the final state under consideration with 
Z — > /i/i and W — > tu t . 

5. ZH : This subprocess can also contribute to the final state Tj + 2 \i (charged— track) 
where the Higgs decaying into a pair of r's, with only one of the r being identified has 
been considered. 

Our chosen event selection criteria have been prompted by all the above backgrounds. 
First of all, we have subjected the events to the following basic cuts: 

• p l ^ track > 10 GeV 

• p^dest-jet > ^ 

• yf**-** 8 > 30 GeV 

• 40 GeV <fi T < 180 GeV 

• \rj\ < 2.5 for Leptons, Jets & Stau 

• AR U > 0.2, ARij > 0.4, where AR 2 = Ar] 2 + A0 2 

• ARfi > 0.2, AR fj > 0.4 

• AR- > 0.7 
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Though the above cuts largely establish the bona fide of a signal event, the background events 
are too numerous to be effectively suppressed by them. One therefore has to use the fact 
that the jets and stau-tracks are all arising from the decays of substantially heavy sparticles. 
This endows them with added degrees of hardness, as compound to jets and muons produced 
in SM process. Thus we can impose a px cut on each track on the muon chamber, and also 
demand a large value of the scalar sum of transverse momenta of all the visible final state 
particles: 

^ pmuonlike track ^ \QQ GcV 

• E\ph\ > 1 TeV 

The justification of these cuts can also be seen from Figures 3 and 4. It may be noted 
that no invariant mass cut on the pair of charged tracks has been imposed. While such a 
cut, too, can suppress the dimuon background, we find it more rewarding to use the scalar 
sum of pt cut. 



3.4 SUSY backgrounds 

Apart from the SM backgrounds, SUSY processes in this scenario itself contribute to the 
final state Tj + 2f + fix + X, which are often more serious than the SM backgrounds. 
These events will survive the kinematic cuts listed in the previous subsection, since they, 
too, originate in heavy sparticles produced in the initial hard scattering. The dominant 
contributions of this kind come from: 

1- X°iX°j production in cascade decay of squarks/gluinos: This is one of the potentially 
dangerous background where both the xT s (hJ = 1? 2) decay into fr-pairs, with only 
one tau being identified. This then mimics our final state in all details with a much 
higher event rate. 

2. u TL Xi production in cascade decay of squarks/gluinos: The decay of v TL as part of the 
cascade produces a Wfi-pair, while the x° decays into a f r-pair to give rise to same final 
state with an additional W which then can decay hadronically. The u t l is produced in 
association with a tau in large number of events (e.g., xt ~~ * ^tlt)- The v TL can also 
be produced from, say, a x\- m both of the above situations, a tau-stau pair can be 
seen together with another stau track, thus leading to a background event. 
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The first background can be reduced partially by looking at the invariant mass dis- 
tribution of the f (having same sign as that of the identified r in the final state) with 
each jet in the final state. If this distribution for any particular combinations falls within 

m o 

m Xi — 20 < Mfj < m Xi + 20 (where m Xi = m x o or — I 2 -, depending on whether x\ or X° * s 
better constructed), we have thrown away that event. The reason for this lies in the obser- 
vations depicted in Figure 2; the invariant mass of a r-induced jet and the f which shows a 
peak close to the mass of the neutralino from which the f-r pair is produced. This has been 
denoted by Cut X in Tables 2 and 3. In addition, if the available information on the effective 
mass tells us that the x% * s better reconstructed in the region, and is produced along with 
a Xi with substantial rate, then a similar invariant mass cut around the X2 m ass will also 
be useful. A further cut on the transverse mass distribution M? Vt (> 1.5m x o or 0.75m x o) 
substantially decreases this background without seriously effecting the signal. 

The background of the second kind can in principle be reduced by vetoing events with 
additional Ws. To identify events with W we have considered only the hadronic decays of 
W's. We first observe the AR separation between the stau (produced in decay of v t l) and 
the direction formed out of the vector sum of the momenta of the two jets produced in W 
decay. If this separation lies within AR = 0.8 and the invariant mass of the two jets lies 
within M w — 20 < Mjj < M w + 20 we discard that particular event. In addition, for a 
sufficiently boosted W , one can have a situation where the two jets merge to form a single 
jet. For such a case, we again look at the stau and each jet within AR < 0.8 around it. The 
invariant mass of the resultant jet is taken to be 20% of the jet energy. The event is rejected 
if a jet with the mass lies within ±20 GeV of the W-mass. We have denoted this by Cut Y 
in Tables 2 and 3. Of course, while it is useful in reducing the background, a fraction of the 
signal events also gets discarded in the process. 

4 Numerical results 

We finally present the numerical results of our study, after imposing the various cuts for all 
the benchmark points. From Tables 2 and 3 one can see that, after demanding a minimum 
hardness of the charged track (p t r f ack > 100 GeV), together with the cut on the scalar sum 
of pt (S|pt| > ITeV), the contribution from the SM processes get reduced substantially. 
The cuts X and Y, defined in the previous subsection, are relatively inconsequential for SM 
processes. However, in the process of solving for neutrino momentum in tau decay, most 
of the SM background events gets eliminated on demanding the invariant mass of the tau, 
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paired with a oppositely charged track, to be around the neutralino mass ijn x o or m x o). This 
is due to the demand that the solution be physical, i.e. the fraction x lies between and 1. 
It is very unlikely to have admissible solutions for x in SM processes, with the r-(muon)track 
pair invariant mass peaking at m x o / m x o . Thus, although the demand < x < 1 is not meant 
specifically for background elimination, it is nonetheless helpful in reducing backgrounds. We 
have verified that the SM contributions within a bin of ±20GeV around the reconstructed 
peak is very small. 

As has been already mentioned, SUSY backgrounds within the model itself is hard to get 
rid of completely. The peak in the transverse mass distribution of the f-v T pair get smeared 
due to such background events (see Figure 4). We have already mentioned two suggested 
cuts, namely, X and Y, which partially reduce these backgrounds. Of these, cut Y suppresses 
(by about 15%) some of the v TL ~Xi/2 events, as can be seen from Tables 2 and 3. the effects 
of this cut on the other SUSY backgrounds as well as the signal are very similar. 

Cut X, meant to eliminate mainly the XrX*j background. Our analysis shows that this 
cut is rather effective in in this respect; the event rate is reduced by almost 50% Surprisingly, 
it also reduces the v TL -Xy 2 background by a considerable amount. The reason for this is the 
following: v TL -x\ii * s produced in cascade decays of squarks and gluinos and the u TL is often 
produced from a xf (the branching fraction being 30% or more in some BP's). In that case 
the decay process is xf ~^ T± v TL - The tau out of such a xt is sometimes identified, whereas 
the tau out of a x?/ 2 (x?/2 — * t± t~ t ) from the other decay chain goes untagged. The invariant 
mass distribution of a track and the jet coming from an unidentified r is clustered around 
m x o (i = 1,2). Thus cut X turns out to be effective in eliminating this type of background. 

After all this effort, however, one still left with background events which smear the peak 
in the transverse mass distribution of the fz/ r -pair. We have to impose an aadditional cut 
on the transverse mass distribution to separate the peak from the background event. This is 
in the form of the demand Mj^ > ^m x o, whereby it is possible to reduce these backgrounds 
further, as can be seen from Tables 2 and 3. It is then possible to determine the chargino 
mass (m> x ±) by looking at the peak, followed by a sharp descent, in the transverse mass 
distribution for several benchmark points. 

The transverse mass distributions for different benchmark points are shown in Figures 5 
and 6. The tau-identification efficiency is assumed to be 50% in Figure 5; Figure 6 reflects 
the improvement achieved in a relatively optimistic situation when this efficiency is 70%. 
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Figure 5: The transverse mass fv T -pair from chargino deacy described in the text, for all 
benchmark points with tau identification efficiency (e T ) = 50%. 
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Figure 6: Same as in Figure 5, but with tau identification efficiency (e r ) = 70%. 
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At BP1 the statistics is very poor and we have realtively few events within a bin of 40 
GeV around m ±. One has about 50% of the events coming from other SUSY processes 
(Table 2). Also the peak is not clearly visible due to the presence of a large number of XiXi 
events, even after imposing the MJ Vt > \m x o cut. 

The situation is similar for BP2 as well. In BP3 and BP4 the peak followed by a sharp 
fall is considerably more distinct, from which one can extract the value of m x ±. 

For BP5 and BP6 the contamination due to the SUSY background is found to be small 
compared to the other benchmark points. 

The xt production rate in cascade decays of squarks/gluinos is also higher there. Hence 
the transverse mass distribution shows a distinct peak, from which a faithful reconstruction 
of chargino mass is possible. 

As a comparison between Figures 5 and 6 shows, the prospect can be improved noticeably 
if one has a better tau identification efficiency (e T = 70%). In such a case, the background 
from XrXi/ X1-X2/X2-X2 i s l ess severe compared to the case where the tau identification 
efficiency is 50%. 

From Figures 5 and 6, one can also see some small peaks in the Mj u distribution with 
very few event rate, in the region where M? Vt > m x o. This can be attributed to those events 
where a X3 or a xl decays into a fr pair, and also to the production of the heavier chargino. 

5 Summary and conclusions 

We have considered a SUSY scenario where the LSP is dominated by a right-sneutrino state, 
while a dominantly right-chiral stau is the NLSP. The stau, being stable on the length scale of 
collider detectors, gives rise to charged tracks, the essence of SUSY signal in such a scenario. 
It is also shown that such a spectrum follows naturaly from a high-scale scenario of universal 
scalar and gaugino masses. 

We have extended our earlier study on the mass reconstruction of non-strongly in- 
teracting superparticles in such cases, by considering final states resulting from the de- 
cays of a xtx\(2) P a i r i n SUSY cascades. The final state under consideration is Tj + 
2f (opposite — sign charged tracks) +^t+X. We have systematically developed a procedure 
for identifying the contribution to pT from the neutrino produced in xt -decay, together with 
a quasi-stable stau. Once this is possible, the transverse mass distribution of the correspond- 
ing f — v T pair can be extracted from data at the LHC, and a sharp edge in that distribution 
yields information on the chargino mass. While eliminating the SM backgrounds in this 
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process is straightforward, we have suggested ways of minimising the contamination of the 
relevant final state from competing processes in the same SUSY scenario. Selecting a number 
of benchmark points in the parameter space, we show in which regions the above procedure 
works. In cases where it does not, the main causes of failure are identified as the overwhelm- 
ingly large contribution from xj-pairs, arid, for example, in the first two benchmark points, 
somewhat poor statistics. The other important issue is the differentiation between the Xi 
and the x° produced in association with the xt- For this, we make use of the assumption of 
gaugino universality as well as the information extacted from the effective mass distribution 
in SUSY processes. 

To conclude, the existence of quasi-stable charged particles, a possibility not too far- 
fetched, opens a new vista in the reconstruction of superparticle masses. We have repeat- 
edly suggested utilisation of this facility in our works on gluino[6] and neutralino[13] mass 
reconstruction. The present work underscores a relatively arduous task in this respect, in 
obtaining transverse mass edges in chargino decays. In spite of rather challenging obstacles 
from underlying SUSY processes, we demonstrate the feasibility of our procedure, which is 
likely to be enhanced by improvement in, for example, the W-and tau-identification efficien- 
cies. 
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Table 2: Number of signal and background events for the Tj + 2f (charged-track)+f^T + X final 
state, considering all SUSY processes, for BP1, BP2 and BPS at an integrated luminosity 300 fb^ 1 
assuming tau identification efficiency e T = 50%. 
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Table 3: Same as in Table-2, but for BP4, BP5 and BP6. 
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